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ABSTRACT 

The same complex matrix model calculates both tachyon scattering for the c = 1 non-critical string 
at the self-dual radius and certain correlation functions of half-BPS operators in N = 4 super- 
Yang-Mills. It is dual to another complex matrix model where the couplings of the first model 
are encoded in the Kontsevich-like variables of the second. The duality between the theories is 
mirrored by the duality of their Feynman diagrams. Analogously to the Hermitian Kontsevich- 
Penner model, the correlation functions of the second model can be written as sums over discrete 
points in subspaces of the moduli space of punctured Riemann surfaces. 
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1 Introduction 



In recent examples of gauge-gravity duality, Feynman graphs of the gauge theory are lifted to 
open string diagrams whose worldsheet holes are summed over, replacing them with closed string 
insertions, to give a closed string theory in a different background. Open-closed string dualities like 
this, such as the 3d Chern-Simons to conifold duality [U [2] or the Kontsevich matrix model to 2d 
topological gravity duality [3J0], have been categorised by Gopakumar as of l F type' [5] because it 
is the faces of Feynman graphs which are replaced by closed string insertions. On the other hand 
dualities such as that between 4d M = 4 super Yang-Mills and Type IIB closed string theory on 
AdS*, x5 5 H are of 'V type' [5] because it is the vertices of the Feynman diagrams, corresponding 
to local operators and interaction vertices, which are replaced by closed string insertions. All known 
open-closed dualities are either of F type or V type [S]. The possibility was raised by Gopakumar 
[S] that for every gauge-gravity duality the closed string theory has open string duals of both types, 
related to each other by graph duality. Topological gravity in 2d was given as an example, with 
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the theory of V type being the double-scaled Hermitian matrix model [3 [9] and that of F type 
being the Kontsevich Hermitian matrix model [3j. Using the proof of equivalence in |10j . the graph 
duality can be shown dynamically by integrating in and out different fields, so that at different 
steps vertices are replaced by faces and vice versa [3]. 

In this paper we show this open-open duality between a complex matrix model of V type called 
the Z model and another complex matrix model of F type called the F model. The Z model is 
known to calculate certain correlation functions of half-BPS operators in J\f = 4 super Yang-Mills 
as well as tachyon scattering amplitudes for the c = 1 string at the self-dual radius. The dual F 
model provides a new Kontsevich-Penner matrix model for these amplitudes and thus throws new 
light on how to write them as integrals over the moduli space of Riemann surfaces. 

The Z model is 



It has two infinite sets of couplings which are often called times in the literature because of the 
relation with the r-function of the Toda integrable hierarchy. 

From the 4d M = 4 super Yang-Mills perspective the Z model is a generating function for 
certain correlation functions of holomorphic and antiholomorphic half-BPS operators built from a 
single complex scalar transforming in the adjoint of the gauge group U(N) I12j. In 'extremal' 
correlation functions, for which the antiholomorphic operators are all at the same spacetime posi- 
tion, the spacetime dependence of the correlation function factors out of the result; the Z model 
computes the remaining combinatorial factor, which is an expansion in 1/JV0 Because local oper- 
ators in N = 4 super Yang-Mills (vertices in the Z model) map to string (or supergravity in this 
case) states, the Z model is of V type. 

In the guise of a normal matrix model |13| the Z model is also a generating function for 
the correlation functions of integer-momentum massless tachyons in the c = 1 non-critical string 
compactified at the self-dual radius. The cosmological constant jjl of the c = 1 string, which controls 
the genus expansion, is related to the rank N of the complex matrix by N = In contrast to 
the double-scaled matrix quantum mechanics (MQM) for the c = 1 stringj^ the Z model requires 
no scaling limit and works at finite N. This is reflected in the fact that the Z model is not a 
triangulation of the Riemann surface itself but rather, through its dual, a triangulation of the 
moduli space of punctured Riemann surfaces. This relation between the half-BPS sector of the 
AdS§ duality and the c = 1 string has been explored in [16^ [T7] based on the similarity of their 
MQM descriptions |18j and shown to be exact at the self-dual radius in [19]. This connection is in 
the spirit of the minimal (p, 1) string embedding in the AdS^ duality |20j . 

The map from tachyons T p with integer momentum p to matrix variables is for k > 



1 This combinatorial factor is unchanged if all the holomorphic operators are also taken to the same spacetime 
position, so really the Z model just generates the two-point function of multi-trace half-BPS states. This is like a 
metric on the multi-trace states, cf. the discussion in [2J • Note too that the extremal correlation functions are known 
not to renormalise when the coupling is non-trivial. 

2 See [141 115] for reviews and references therein. 
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The individual tachyon correlation functions are then 



Tk x --- T kp 71^ ■T 



9/ C=l 



tr(Z fcl ) • • • tr(Z k ?) tr(Z tfcl ) • • • tr(Z tfc ') 



(3) 



On the righthand side the correlation function is taken using the complex matrix model with 
Gaussian action tr(ZZt). It is computed by Wick-contracting with the propagator 



Z) zt») = 



(4) 



These correlation functions can be computed to all orders in N using symmetric group techniques 

021 H]. 

The dual F model is exactly the same function of {t}, {t} 



Z({t},{t}) =J[dF]^ xn e -*(W)-Ntrlo g (l-A-iFB-W) 



•tr(FFt)+JV E^Li ^tr[(A- 1 FB~ 1 Ft) fe 



(5) 



The couplings {t} and {t} are encoded in matrices A, B by a Kontsevich-Miwa transformation 

i =V— 1 



i=l 



kaf? k 



trA~ 



tk ^ kb k - k UB 

3=1 3 



(6) 



To expand Z({t},{t}) in these variables, compute correlation functions with the even-valency 
vertices that appear in ([5]) for k > 1 using the propagator from the matrix model 



p 1 * ptfc 
3 1 



1 



(7) 



The colour index for each face of the F model Feynman diagrams comes with either an Oj or a bj , 
so the couplings {t} and {t} are associated to faces of the F model. Thus the F model is of F type. 

Relations between Hermitian matrix models via graph duality have appeared before in the 
literature, as have complex matrix models similar to the F model (see for example }23]-|32j. also 
in connection with r-functions \33\ 134]). The F model is of Kontsevich type because the couplings 
are encoded and expanded similarly to the Kontsevich model for topological gravity [3]. It is also 
of Penner type because the appearance of a logarithmic term in the action is similar to the Penner 
model for the virtual Euler characteristic of the moduli space of punctured Riemann surfaces M g , n 
[35j . In fact the F model is a complex matrix model analogue of the Hermitian Kontsevich-Penner 
model studied by Chekhov and Makeenko [36], which is dual to the Hermitian version of the Z 
model (before the double-scaling limit) in exactly the same way [55| 15]. 

The most direct way to prove the duality between the Z and F models is using character 
expansions, see Section 21 Term-by-term it can be seen that the Feynman diagrams of the different 
models are graph-dual. In Section 14.41 it is shown that the correlation functions are in fact sums 
over Hurwitz numbers, which count holomorphic maps from the worldsheet to CP 1 branched just 
three times. The valencies of the vertices and faces specify the ramification profiles. 

Another proof is given in Section [2] with the techniques used in the 2d topological gravity case 
by Maldacena, Moore, Seiberg, Shih [TU] and Gopakumar [5], which involve integrating fields in and 
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out twice. In this method the graph duality between the Z and F model can be seen 'dynamically', 
as explained in Section [3j Every Feynman diagram in the original Z matrix model corresponds to 
a diagram in the F matrix model to which it is dual. This insight is crucial to read off the correct 
terms that are identified in the different models. In fact with the propagator the F model is 
only sensitive to 'skeleton' graphs of the Z model where propagators running parallel between the 
same vertices are bunched together into the same edge. These skeleton graphs were introduced 
in [38] as part of Gopakumar's programme to find the closed string duals of free gauge theories 
[371 ESI ES]. 

An advantage of the F model is that its correlation functions can be expressed directly as 
integrals over the moduli space Mg,n of punctured Riemann surfaces, using the example set by 
the Kontsevich model [3]. In the Schwinger parameterisation of the propagators, the Schwinger 
lengths associated to each edge of each Feynman graph provide coordinates on a cell decomposition 
of M g ,n- The integrals over the top-dimensional cells in M g , n require all vertices of the graphs to 
be trivalent. The vertices of the F model have a minimum valency of four, which means that the 
correlation functions can only come from lower-dimensional cells in the moduli space. Furthermore, 
following the analysis of the Hermitian Kontsevich- Penner model in [30], the integral localises on 
discrete points in these subspaces, see Section [5j 

Despite the fact that the Z model needs no double-scaling limit for its identification with the 
c = 1 string, it is still possible to take one. From the N = 4 perspective it is the BMN limit [UJ and 
it limits the F model to only 4-valent vertices, see Section 14.31 The meaning of this limit for the 
c = 1 string is unclear. Rewriting the correlation functions of the Z/F model in terms of Hurwitz 
numbers in Section 14.41 this limit involves restricting to a special class of Hurwitz numbers called 
double Hurwitz numbers, with arbitrary branching profiles at two points and simple branchings 
elsewhere. 

Another topological matrix model for the c = 1,R = 1 string is the model of [421 143] . 
reviewed in [44] . where just the positive momentum tachyon couplings are rearranged in this way 

2({t},{t}) = J [dM] H NxN e M-M+E£^(^-)*) (8) 

The integral is over a Hermitian matrix M. The relation of the Z model to this model was 
explained by Mukherjee and Mukhi in |45j : a direct transformation of the model into the F 
model is shown in Appendix O 

2 Proof of duality using integration in-out-in-out 

In this section the duality between the Z and F models is proved using the techniques of [101 [5] by 
integrating in and out different fields. This makes the graph duality of the models manifest, as is 
explained in the next section. 

The partition function for the Z model is 

Z({t},{t}) = J [dZf NxN e-M^+ESLifcMzfrj+E^M^*) ( 9 ) 

This model is the same as the Model II for the c = 1 string at the self-dual radius R = 1 with 
N = v = — ifi in [13] • Although the integration in [13J is over a normal matrix with the condition 
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[Z, Z*] = enforced, with this action for R = 1 both the complex and normal matrix model are the 
sameo The tachyon scattering matrix agrees with older results calculated in the literature |47| , 

Substitute the tk and if. for two diagonal n x n matrices A and B, with eigenvalues and bj 
respectively, using the Kontsevich-Miwa transformation 



1 

k 



— k 



1 



n 



n 



1 



(10) 



For the tk to be independent whenever the ti(Z k ) are, we need n > N and similarly for the tk- 

The exponentiated tr(Z k ) operators can be written as inverse determinants provided the a, are 
sufficiently large (to avoid convergence issues) 



exp 



X>tr(Z A 



k=l 



exp 



n . 



,fe=l i=l 



exp 



^ tr log I 1 



i=i 



Z 



n 

8=1 



det 1 



(11) 



[In the 2d topological gravity case the determinants in the double-scaled Hermitian matrix model 
correspond to exponentiated macroscopic loop operators for FZZT branes 



trlog(aj — M) = J ^-tre~ 



l(a,-M) 



(12) 



Each of the n FZZT branes has boundary cosmological constant a^. There is no clear such inter- 
pretation of the determinants as wavefunctions of FZZT branes here, and in fact the more natural 
extension of [10] would be to investigate macroscopic loops in the matrix quantum mechanics, cf. 
[49] . The fact that we have inverse determinants in the c = 1 case (also present in the study of the 
normal matrix model in [35]) also differs from [10] and alters the statistics for the fields that we 
integrate in later, which in the c < 1 case are fermionic strings stretching between the ZZ [50] and 
FZZT branesQ From the 4d M = 4 SYM perspective, these determinants (more clearly expanded 
in equation (|48p ) are interpreted as giant graviton branes in the bulk |5H 118] .] 
Using ([lip the Z model partition function is now 



Z{{t},{t})= [dZ] 



INxN 



- tr(ZZt)+Er=i h tr(Z fe )+£r=i % tr(Zt*) 



det ( ^det (B) ^n 5i ^ 



Zt) 



(13) 



The correlation function is taken with the Gaussian action tr(ZZ'), as will always be the case for 
the Z model. 



3 A normal matrix can be decomposed into a unitary matrix U and a diagonal matrix of its complex eigenvalues 
D, Zn = UDU^ . This is not true for a complex matrix, for which we have Z = U(D + R)U^ where R is strictly upper 
triangular [46]. It can be checked that in the action (j9j R completely decouples, and since the measure on U and D 
is the same, the normal matrix model is equivalent to the complex matrix model with this action. See equation (|92p 
in Appendix [C] for an alternative way to decompose a complex matrix. 

4 Note that if we had chosen to include a minus sign in the identification (|10|l we would have had normal determi- 
nants here and fermions integrated in later. The choice of sign is left to a physical interpretation in the future. 
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Writing the products of determinants using single determinants of larger nN x nN matrices, 
we can write them as integrals over two sets of complex bosonic fields^ 

1 1 



det(A ® % - In, <g> Z) det(B ® I N - l n ® Zt) 

[dZ] NxN [dC] Nxn [dD] Nxn e -tr[ZZ^CHAm N -I^Z)C + DHB& N -l n ^)D] (M) 



C e i and D e j are bifundamental fields with e = 1, . . . N and i,j = 1, . . . n. Again because we have 
inverse determinants this contrasts to the minimal string case, in which one must integrate in 
fermions rather than bosons. 

Next integrate out the Z field, after rewriting (|14p appropriately 

J ' [dZ][dC][dD] e -*[&-DD*)(z i -C&)-CCftDDl+&AC+DiBD] 

= J[dC][dD] e -tr(ctAC+Dt B D-CCtDDt) (15) 

This is the C, D matrix model. The quartic vertex is CC^DD^ = C e iC\-D fjDl-. It has propagators 

The C, 15 model can also be expanded as a function of the couplings {t}, {t}, see Appendix [Bl 
Next we integrate back in an n x n complex matrix F % - being careful with the indices 

J [dC][dD] e ~ C ei 1 A h i 2 CEi 2~ D lj 1 B ilj2 D ^2 +CEiD lj D f3 C }i 

= J [dF]nxn[dC] [dD] e -t F )- c ^ D l^ F ^- D ii c )^ c ^ 

= J [dFf nxn [dC][dD] e -tr(FFt-D^tC-C^D + C^AC + D^BD) (17) 

To integrate out C and D write them together as a single N x 2n field so that the cubic terms 
become 

The result is an inverse determinant of an [N x 2n] x [N x 2n] matrix 



-1 f / \ \ —AT 

4 —F \ I 



/ [dFl jdet ( / ) ® I N I e- tr ^ + ) = | [dF] jdet ( 

Now write the determinant in terms of a product of matrices 



-Ft B 



•tr(FFt) 



dCt -F» "/)= det fl H-bV "T' 



d e t(/l)det(B)det( _ B 1 1J?t /4 1 ' F ) < 19 > 



' This type of identity for the determinants only works with the determinant of a Hermitian matrix. We extend 
it to our case by noting that the Z model only depends on the eigenvalues of Z and since they are complex we can 
extend the identity by analytic continuation. 
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The constant terms det(A) det(B) cancel those in (fT3j) . so we get 



Z({t},{t})= [dF] e 



■tr(FFt)-ATtrl 



1 



-A~ X F 
1 



(20) 



Expanding the logarithm for large A, B, just as we did in (jlip . we find 



Z({th{t})= [dFt xn e 



-tr(FFt)+ATJ]^ =1 Itr 



A~ l F 
B- X F^ 



J [dF]^ xn e - tT ( FF ^ +N Zk=it tT [( A ~ lFB ~ lFt )' 



(21) 



This is the F model introduced in equation ([5]). To extract the propagator study the quadratic 
term 

F] to - NiA-^KB- 1 )^ Ftf (22) 



Using (A 



a a 



8\ the propagator is 



3 1 



(1 - Na^bJ 1 ) 

Alternatively we could have taken the plain quadratic term tr(FF') with plain propagator 



(23) 



FjF**) 



I pla 



(24) 



and treated the k = 1 term TV tT{A- l FB~ l F^) from §FQ as an additional interaction vertex. The 
propagator (|23p is then a sum over an arbitrary number of intervening such 2-valent vertices 



NP 



6j6* 



— - ( F] tv{A~ l FB~ l F^) fV) =Y(Na7 1 bJ 1 )P Sis'? = - 
^ «' \ 3 /plain ^ 1 1 ' 1 3 (l-Na^b' 1 ) 



(25) 



p=0 r p=0 

This is an important issue for the interpretation of the graph duality in Section [3J since these 
2-valent vertices are exactly those which are dual to the faces bounded by parallel propagators 
between the same vertices, see Figure [7J Bunching such parallel propagators into a single edge, 
reducing the Feynman diagram to a skeleton graph, removes such 2-valent vertices from the dual 
graph. This works because the propagator (|23p sums over all possible Z diagrams with the same 
skeleton graph. These interpretational issues are also crucial to understand the dual of the planar 
two-point function of the Z model and the character expansion to which we turn in Section [H 
It will often be useful for calculations to rescale the F model F — > \J~AF\J~B to get 

Z({t},{t}) = (detAHdetB)" J [dF]^ xn e-^^+^^iH^f (26) 

This model then has propagator 



F] Ft*\ 

/ scaled OjO 



N 



(27) 



This form is useful for transferring the Hermitian Kontsevich-Penner analysis of [50] to the F model 
in Section [SJ Example F model correlation functions are computed in Appendix [XI 



S 



3 Dynamical graph duality 



[All figures are printed at the end of the paper to avoid cluttering the text.] 

In this section the duality between the Z model and the F model is shown to work at the level 
of individual Feynman diagrams. An explicit proof with formulae is provided in Section HJ 

Each graph of the Z model corresponds to a graph of the F model which is dual to the original 
Z model graph. The graph duality is shown 'dynamically' in the sense that it is split up into stages 
where we first replace the vertices of the Z model by faces (integrating in C,D), then contract and 
expand propagators in different channels (integrating out Z and in F) and then finally replacing the 
faces of the Z model by vertices of the F model (integrating out C,D). This analysis completely 
mirrors the 2d topological gravity analysis by Gopakumar [5]. 

We start with a correlation function such as ([3|) for the Z matrix model with several holomorphic 
vertices tv(Z k ) and several antiholomorphic vertices ti{Z^ k ). Each possible Wick contraction with 
the propagator will give a different Feynman graph. Because the only non-trivial propagator is 
Z — )• Z' and no vertices mix Z's with Z''s, there can be no propagators connecting a vertex back 
to itself. The Feynman diagram has a minimum genus surface on which it can be drawn with 
no lines crossing. Propagators that run parallel to each other between the same pair of vertices 
will be refered to as 'homotopic' and can be bunched together into a single edge. If we do this 
for all possible propagators then we get what Gopakumar named a skeleton graph |38j . For each 
correlation function there may be several topologically distinct skeleton graphs, each corresponding 
to a number of possible Wick contractions of the original correlation function. Later in Section H] 
we will see that these distinct skeleton graphs correspond to the cut-and-join operators of [21] . 

It should also be noted that some graphs, where there is more than one single trace for both 
holomorphic and antiholomorphic operators, are disconnected; cf. the examples in Appendix I A.4i 

When we integrate in the C, D matrices, a vertex like tr(Z m ) is replaced by a face of C's and 
an antiholomorphic vertex tr(Zt m ) is replaced by a face of D's, see Figure [TJ In the figures the 
double lines of the Z are drawn with solid lines e = 1, . . . N, while the bifundamental C's involve 
both a solid line and a dashed line for i = 1, . . . n and the D's a solid line and a dotted line for 
j = 1, . . . n. We get cubic couplings C^ZC and L>t Z^D from the action in (|14p . 

We then integrate out Z, Z^ . The Z — > Z* propagators shrink to give graphs with the quartic 
vertex CC^DD\ cf. Figure [2j This is the C, D matrix model from (|15p . It can be expanded in its 
own right, see the examples in Appendix [Bl 

Integrating in the complex matrix F, the quartic vertex CC^DD^ expands in a different channel 
into an F — > Ft propagator with the cubic couplings C^FD and D^F^C of (|17p at each end, cf. 
Figure In this way each propagator of the Z model corresponds to a transverse propagator of 
the F model via the quartic CC^DD^ vertex. 

Finally we integrate out C and D to get the dual graph in terms of F, Ft for the F model 
(|2ip . Each face involving a solid line, corresponding to the original faces of the Z model Feynman 
diagram, is replaced by an even- valency F, Ft vertex (FF^) P , cf. Figure HI Faces to which only two 
F ->• Ft propagators connect become just a sing le F Ft propagator. These faces correspond 
to the faces between parallel homotopic propagators of the Z model, cf. Figure UJ Thus a bunch 
of parallel propagators from the Z model become just a single F propagator in the dual graph; in 
other words the F model is sensitive only to the topology of the skeleton graph. 

We have thus seen that the correspondence between the Z model and the F model corresponds 
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to graph duality. A vertex of the Z model becomes a face of the F model; homotopically bunched 
propagators (edges) of the Z model become a single propagator (edge) of the F model which is 
perpendicular to the orig inal Z -> Zt propagators; and finally the faces of the Z model become 
even-valency vertices of the F model. 

An important constraining feature of the F model (|21|) is that the faces of the graphs are always 
associated with either A's or B's but never both. Prom the propagator in (|23p the a% is associated 
to one (dashed) index line while the bj is associated to the other (dotted) index line, cf. the left 
part of Figure [5j The vertices also preserve the Oj and bj associations, cf. the right part of Figure [5j 
This is because the faces correspond to vertices of the Z model where the A's map to tr(Z k ) vertices 
and the B's to tv(Z^ k ) vertices. In fact, because a Z model vertex can never have a self-contraction, 
each edge of the F model has an A face on one side and a B face on the other, reflected in the fact 
that the double-line propagator for F has one dashed and one dotted line. 

A full example of this dynamical graph duality is given for (tx(Z 2 ) ti(Z) tr(Z^ 3 )) in Figure El 
To conclude this section we summarise the procedure for seeing the duality from graphs: 

• Starting from the Z model, expand the partition function in correlation functions of holo- 
morphic and antiholomorphic operators. Each correlation function corresponds to a sum of 
different topologically-distinct skeleton graphs. Each skeleton graph corresponds to a sum of 
topologically-identical F model Feynman diagrams, to which the skeleton graph is dual. 

• Starting from the F model, expand the partition function in correlation functions of the 
interaction vertices. Each F correlation function splits into classes of topologically-identical 
Feynman diagrams. Each such class of Feynman diagrams is dual to a family of topologically- 
identical Z skeleton graphs. 

We will now give an alternative proof of this graph duality using symmetric group techniques. 



4 Character expansions 

In this section we use symmetric group and representation theory techniques to expand both the 
Z and F models and show that they are equal. The graph duality is revealed in the character 
expansions of the models. In fact a more general correspondence is proven 



Z({t},{t},{s}) = J[dZf NxN e -MS-^s-'^)+Er=i* fc tr(^)+ESLii fc tr(^) (2g) 

^2 \R\ ] - xr(S)xr(A- 1 ) Xr (b- 1 ) (29) 

l(R)<N ^ R 

= J[dFf nxn e-trFjrt+E^.fctrCX-^B-iFt)* (3Q) 

Ocj[dF]^ xn e-^AFBFfy+^s^F^ (31) 

In the first line the quadratic term of the Z model has been modified with an N x N diagonal 
matrix S whose eigenvalues are S± - ■ ■ SV- The propagator is now 

Z) Z^ h ) = sis] 5 e h 5 9 f (32) 
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Each face of the Z model is bounded by 2k edges (the face has valency 2k). Since each edge picks 
up a factor of S? from the propagator, each 2A:-valent face of the Z model comes with a factor 
tr(S k ). Each vertex of the Z model has a coupling tk or tk depending on whether it is holomorphic 
or antiholomorphic. To recover the vanilla Z model of ([JJ set S to be the identity matrix S = Ijy. 

In the second line the Z model has been expanded in characters xr °f the group U(N) whose 
representations R are labelled by Young diagrams with at most N rows, 1{R) < N. \R\ is the number 
of boxes in the Young diagram and dn is the dimension of the symmetric group representation for 
the same Young diagram. If we set £ = Ijv here, the character of S becomes the dimension of the 
U{N) representation R, xr(Jn) = dimjv-R- Note the symmetry of the character expansion in the 
variables {£}, {t} and {s}, as encoded in the matrices A~ x , B~ l and S. 

In the third line this character expansion is equal to the F model, which has been modified 
so that now each 2/c-valent vertex is weighted with a coupling S}.. This coupling is related to the 
N x N matrix S by the Kontsevich-Miwa transformation s& = r tr(S k ). If we set S = In here, Sk 
becomes Sfc = r tr(I^) = and we recover the vanilla F model of (JSJ. 

In the fourth line the F model has been rescaled F — > \J~AF\f~B as in equation (|26p to make 
the duality with the Z model in (|28p more transparent. 

In the Z model each 2/c-valent face picks up a coupling Sk from the propagator while the fc-valent 
vertices pick up couplings tf~ or tk- On the other hand, for the dual F model, exactly the reverse 
happens: it is the 2/c-valent vertices which pick up the couplings Sk while the fc-valent faces pick 
up couplings tk or tk from the propagator. Thus the addition of the matrix S and its associated 
couplings Sk allow us to follow the graph duality more clearly. 

Similar Hermitian models where faces of their graph expansions are weighted separately from 
the vertices have been studied in [U]-[27] and indeed the F model and its character expansion were 
mentioned in the review [25] as a complex generalisation. 

4.1 Character expansion of Z model 

In this section the symmetric group techniques of |18} |2"T] are used to expand the Z model, which 
generates the correlation functions computed in those papers. 

The first step is to write each product of holomorphic traces (often referred to as multi-trace 
operators) using a conjugacy class element of Sk where k is the sum of the powers of Z . We write 
the conjugacy class as a partition /X2, • • • n p ] of k so that k = J2i Mi = Each m corresponds 
to a ^j-cycle in Sk- In this way 

tr(Z^) tr(Z^) • • • tr(Z^) = Z^Z^ ■ ■ • Z^ = tr(a Z® k ) (33) 

where a is in the conjugacy class [fj,\, fj,2, • ■ • Mp]- For a partition [i of k (written fi h k) we write a 
representative of the corresponding conjugacy class \p] as a M G \p] C Sk- In fact for these operators 
it does not matter which representative we pick because all elements in the conjugacy class [//], 
given by p~ l a^p for p G Sk, correspond to the same multi-trace operator (f33|) . The size of the 
conjugacy class is 

M= Ml = Y^~^HM\ (34) 

6 Assume the parts of the partition are ordered fit > fii+i so that the parts map to the row lengths of a Young 
diagram with k boxes. 
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|Sym([/z])| is the size of the symmetry group Sym([/i]) of the conjugacy class. Sym(a At ) is the 
subgroup of Sk that leaves a M invariant under conjugation: Sym(a M ) = {p S S). \ pa^p" 1 = a M }. 
i p (p) is the number of parts of p of length p. The factors in the denominator of (|34p come from 
the cyclic symmetry = Z p of the i p (p) cycles of length p and the i p (p)\ ways of exchanging them. 

Summing over all permutations of Wick contractions of the fields, the correlation functions of 
the vanilla Z model ([I]) are an expansion in N where the power of N is the number of faces of the 
Z model graph [18] 

(tr(a v Z^ k ) tr(cv Z® k )^ = ^ N C(a,ra^) = N^diaa^ra^T- 1 ) (35) 

C(a) is the number of cycles in the permutation a £ Sk (i.e. the number of parts of the corre- 
sponding partition). The second equality is just a rewriting of the first with the symmetric group 
delta function <5(c). 5(a) is zero on all elements of Sk except the group identity id, on which it is 1. 

Since the operators only depend on the conjugacy of it is often useful to write them with 
sums over the entire conjugacy class tr(a^ Z® k ) = tAt tr(£^ Z® k ) where 

s M = y, a ( 36 ) 

ae[M]c5 fc 

The sum over a £ Sk in ([35]) can also be sub-divided into conjugacy classes 

J2N C ^ a = Y,N Ci[X]) Zx (37) 

aeSk Ahfc 

This element of the group algebra is often written in the literature as 

^JV-*£^ (W)Ej = E ^_ Sa (38) 

Ahfc Ahfc 

We have used the identity T([A]) = k—C([X\), where T([A]) is the minimum number of transpositions 
needed to build an element in the conjugacy class [A]. 
The correlation function ()35[) can now be written 

(trK Z^ k ) tr(a, Z® k )) = *(E„n*E M ) = pjpj Y N ° m S^WJ (39) 

Each summand in (|39|) now corresponds to topologically different graphs. The partitions p, v label 
the different holomorphic and antiholomorphic vertices respectively. The partition A, on the other 
hand, labels the faces. Each part \ of the partition labels a face so that C([A]) is the total number 
of faces and hence gives the power of N. Putting a non-trivial matrix S into the Z model as in 
pgj) the factor N C W> is refined 

k 

N c([\]) ^ tr(aA g»k) = JJ [tr^jjipW ( 40 ) 

v=\ 

i p (\) is the number of parts in A of length p. Each part of length p corresponds to a 2p-valent 
face of the Z graph bounded by 2p edges. The parts of length 1 correspond to faces with only 
two edges; these are the faces between propagators running parallel between two vertices, see the 
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lefthand diagram in Figure [7J Bunching these parallel propagators into single edges so that we 
get a skeleton graph corresponds to ignoring the parts of A of length 1. This proves a conjecture 
in [21] that the different E A correspond to the different skeleton graphs if you ignore the parts of 
length 1. The genus of the Feynman diagram can be read off using the Euler characteristic with 
V = C([fi]) + C{[v]) vertices, E = k edges and F = C([A]) faces 

x = 2 -2g = V-E + F = CM) + C{[u]) + C([A]) - k (41) 

In |21] the E A were called cut-and-join operators because they have an action on a multi-trace 
operator tr(S At Z® k ) that can split single traces into many traces or join many into one. This action 
is simply left-multiplication 

E A : tr(S M Z® k ) i ^ tr(E A E M Z® k ) (42) 

For example the conjugacy class of transpositions A = [2, l fc ~ 2 ] gives the cut-and-join operator E™. 
Acting on the single trace operator given by /i = [k] splits it into all possible double-traces \p, k — p] 

E [2] tr(Z fc ) = ^tr(E M E [Jb] Z® k ) = ^te(Z*)ta(Z*-*) (43) 

Acting on a double trace with E[ 2 ] can also join it into a single trace. Similar results follow 
when acting with more general cut-and-join operators E A on more general trace structures. When 
writing partitions A for cut-and-join operators we omit parts of length 1 so that in this example 
A = [2, l k ~ 2 ] — > [2]. This omission of length- 1 parts corresponds to only considering skeleton graphs, 
as discussed aboveQ 

Next we turn to the character expansion. The symmetric group delta function in ([39]) can 
be expanded as a sum over representations R of Sj-, also indexed by partitions, in terms of the 
characters XR( a ) °f a m the representations R and the Sk representation dimensions (Ir = x#(id). 
The identity is 5(a) = ^^m-k^RXR^)- Using the fact that d R x R (H^a) = xr{v)xr{^> y) , since 
E^ is central in S^, and the relation between an d the unitary group U(N) dimension of the 
partition R, N k XR{^k) = k\ dimjv -Ru the correlation function (|39|) has a character expansion 

tr(a„ Z^ k ) tr(a„ = V fc! XrMxrM R 

Rrk n 

This identity can also be seen by expanding the traces in terms of U(N) characters tr(a^ Z® k ) = 
Tl<R!rk XR( a y)XR(Z) where the U(N) character xr(%) is itself expanded in terms of traces xr{%) = 
M SaeSfe X-R( a )tr(a Z® k ). The two-point function of these characters xr(Z) 1S diagonal [18]. Note 
that the dimension dimjy R vanishes if the number of parts l(R) of the partition R exceeds N. 
For general S as in fj28[> the correlation functions ()39p and (|44p become 



(tr(a„ Z^ k ) tr(a„ Z® k )) = £ £ |[A]| tr(a A S ^) Xr{ ^ )Xr ^ )xrM (45) 

XR 



E k\ XR{uii)XR{av)XR(S) 
r] r> 



Rrk 



7 In [21] it was useful to think of the correlation function (|39|) as an inner product of tr(a^ Z'® k ) with the result 
of Qk acting on tr(a jU Z® k ). The bra-ket notation of the inner product in [21] is the same as the delta function here 
(u\y) = w ^- Tl S(E^) = |Sym(M)|^=„. 

8 Characters of the cut-and-join operators X-r(^k) roughly correspond to U(N) Casimirs, cf. Section 2.7 of [3T] 
and original references in [53] \Mi ■ 
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Now we are in a position to expand the full Z partition function. It is useful to collect the 
expansions of the exponentials into multi-trace operators indexed by partitions \x of k, so that the 
total number of Z fields in each term is k. Using the Kontsevich-Miwa transformation (|10p 

The denominator factors give the size of the conjugacy class from (|34|) so 



OO i r i i oo 

7feN 



e Er =1 * fc tr(^) = ££^tr(a M (^J**) tr(«, = ^-£xr(A-1) X r(Z) (48) 

k=0 tihk ' k=0Rhk 



The Z model partition function (|28p becomes 

z({t},{t}) = L^w* k ) e Er=i**^») 



= E E TT tr ^ ^" 1 ) 0fc ) ^T tr ^ (-B -1 )^) (tr(^ zt® fc ) tr(a M 

k=0 fj,,uhk 

j2^2 kl xg^Mg^Mg) (49) 

fc=0 i?hfc ^ 

Noting that Xr{S) f° r an. N x N matrix S 1 vanishes if Z(i2) > iV, and since n > N this constraint 
takes precedence over those coming from xr^A" 1 ) and xr,{B~ 1 ), this concludes the proof of ([29|) , 

4.2 Character expansion of F model 

A generalisation of the Itzykson-Zuber integral [22] (see for example [32]) can be used for the 
following integral for the F model with plain quadratic term (|24p 

XniA^FB-Wj) = f [dFf nxn e -^ FF \ R (A^FB^F^ = ^ ^ A ^ B ^ (50 ) 

/ plain J dR 

For an F model correlation function the vertices are dictated by a partition A of k, so that the i p (X) 
parts of A of length p correspond to i p (\) 2p-valent vertices. Expanding with the formula (|50p 

TT |tr(A- 1 FB- 1 F t ) p l ) = ( tr(a A (A~ 1 FB~ 1 F^)® k ) 

~ L J / plain \ / plain 



p=l 



= E ^ (51) 

Rrk n 

Expanding out the characters xr(^4 _1 ) and \r{B~ 1 ) as multi-traces gives 

E E tr(a„ (A" 1 )**) tr(a„ (iT 1 )®*) ^^)x^K)x fi K) (g2) 

ill- A; fi,v\-k 

This expansion is important because we see how the faces of the F model diagrams are controlled 
by the partitions fj, and v that organise the couplings {t} and {£}. Comparing (|52p to its Z model 
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equivalent (|45p and picking out the terms for particular face configurations in the respective models 



k\ (tr(£ A (A^FB^F^ fP)) = (tr(I!„ Z^ k ) tr^ Z® fc ) 

\ / plain, ^t,i^ term \ /A term 



Rhk n 

= kl SpuYixVn) (53) 

In other words both models are doing nothing other than computing the conjugacy class algebra 
of the symmetric group. 

If we use the corrected propagator (|23|) then there are no 2-valent vertices in the F model, 
corresponding to the fact that all parallel propagators of the Z model are bunched into single 
edges. Thus with the corrected propagator for the F model the partitions A of vertices have no 
parts of length 1, i\(X) = 0. To relate this to the above analysis above define A r = A + [l r ] so that 

oo 

k+r 



fc!(tr(£ A (A- 1 FB- 1 F^) 9k )\ = kl (tr(£ A r (A~ l FB~ l F^) 



plain 

r=0 

oo 

= 12 12 tr ( a » (A~ r f k+r ) tr(a„ (B- 1 )®^) (tr(p v Z^ k+r ) tr(S M Z® k+r \ 

r=Q )i,u\-k+r 

Even if we fix the structure of //, v, A there are still an infinite number of Z diagrams for every F 
diagram because in the Z model we can have any number of bunched parallel propagators at each 
edge, corresponding to adding arbitrary numbers of 2-valent vertices to the plain F model. 

In the original vanilla Z model with S = Ijy picking a particular power of N in the expansion 
of the correlation function (|39p corresponds to fixing the genus of the diagrams in the Z model 
expansion. Looking at the expansion of in terms of cut-and-join operators £ A in equation 
(|38p . fixing the genus means we only consider A composed from a fixed number of transpositions 
T([A]) = p. The generic A with this property is [2 P ], followed by cases where the transpositions 
blend into longer cycles [3, 2 P ~ 2 ], [4, 2 P ~ 3 ], [3, 3, 2 P_4 ] etc. For example [2,2] — > [3] corresponds to 
the degenerate multiplication of two transpositions when they share an element i to get a 3-cycle 
(ik)(ij) = (ijk). 

In the F model A = [2 P ] means p 4-valent vertices in the graph. The other terms for fixed genus, 
such as [3,2 P ~ 2 ], correspond to vertices in the F model colliding to create higher valency vertices. 
For example the collision of two four-valent vertices into a single six- valency vertex for the dual to 
the Z model torus two-point function in Figure [9] is the counterpart of [2,2] — > [3]. The limit with 
only the generic cut-and-join operator S[ 2 p], corresponding to having only 4-valent vertices in the 
F model, is discussed in Section 14.31 

Note that the number of edges E in the Z or F model diagram is given by a sum over the 
faces/vertices weighted by their valency 2E = J2p=i^PhW = Yl p =iP (*?(/•*) + hi 11 ))- This is 
automatic since A, fi and v are each partitions of E = k. 

Finally, to prove equation (|3U|) for the full expansion of the F model partition function expand 
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the exponential using ([35]) . With couplings s& = \ tv(S k ) to the 2&-valent vertices 
= f tr(a A S® fe ) / tr(a A (A _1 i ?1 -B _1 F t )® fc ) , 

^— ' ^— ' k\ \ I plain 



k=0 Xhk 



_j2J2 \R\ ] -xr(S)xr(A~ 1 )xr(B- 1 ) (54) 
4.3 A large k limit 

It was shown in |55[ [21] that if you take a given correlation function in the Z model and fix the 
genus, corresponding to a particular power of N, then in the large k limit of operators with many 
powers of Z the cut-and-join operator Xr 2 pi dominates (cf. the torus two-point function in (|75p ) 

ti{a u Z^ k ) tr(«^® fe )) = TW7T E ^ *( S " S a s m) 

P / JV^-P term U M ^ P 

~*° MM (55) 

In other words, if we take the expansion of fi^ in cut-and-join operators in equation ()38p and take 
this limit for each inverse power of N then the £[ 2 p] term always dominates as a function of k in 
the correlation functions at this genus 

k ^ oo 1 

«* = E^ E ^ -~ E^V] (56) 

p=0 A:T(A)=p p=0 

For a fixed genus the terms we lose in this limit correspond to the degenerate cases where transpo- 
sitions collide to give higher cycles. These degenerate cases appear with less frequency so they are 
suppressed for large k. 

Since in this large k limit Spp] ~ pT^pl' ^ k can " 3e sa ^ ^° ex P onen tiate $7^ — > exp(-^S[ 2 ])- 
The sum over transpositions S[ 2 ] has a simple interpretation in terms of either splitting or joining 
a trace. Geometrically this means that n punctured genus g surfaces factorise into 3-punctured 
spheres. See [2T] for more details. 

Here N has been treated as a book-keeping device for the genus. If we allow N to be a 
number, then the limit must be taken delicately. Really we are taking the double-scaling BMN 
limit k ~ — y oo [4JQ and expanding with the non-planar coupling </ 2 = jr < 1, cf. [21]. If 
k grows faster than then other terms become significant and spoil the simple exponentiation 
n k -> exp (^£[2])- 

In terms of the F model, restricting to cut-and-join operators of the form S[ 2 p] means that we 
only have 4- vertices in the model. In terms of the couplings Sf~, we are setting s& = for k > 3 
and only retaining 4- vertices s\ = N, s 2 = y. Diagrams involving higher- valency vertices resulting 
from collisions of 4- vertices are dropped. 

In the 2d topological gravity case the Hermitian matrix model with couplings X}fc=i tk tr(M fc ) 
can also be rearranged exactly into a Kontsevich-Penner model before taking the double-scaling 
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limit. This is the Chekhov-Makeenko model [36] and is the Hermitian version of the complex 
matrix model duality discussed here. The dual model has vertices of all valency, including odd 
ones. Taking the double-scaling limit eliminates all but the 3-valent terms, giving the standard 
Kontsevich model. 

4.4 Hurwitz theory for CP 1 \ {0, oo} 

The formulae for the correlation functions as sums over cut-and-join operators (|39p can be inter- 
preted in terms of Hurwitz numbers that count holomorphic maps from the Riemann surface E> g on 
which a graph is drawn to the sphere CP 1 with three branch points. For a A;-sheeted covering of 
CP 1 by a genus g surface with ramification profiles /i, v and A the number of coverings is 

Cov^,A,/,):= £ = ^(S,E A E M ) (57) 

The maps are counted up to automorphisms of the covering map / : S g — > CP 1 . The genus g for 
which this is non-vanishing is given by the Riemann-Hurwitz theorem, which relates the genus g of 
E) g to the branching numbers at each branch point 

2g - 2 = -2k + T([u]) + T([A]) + T(\jm]) (58) 

The branching number at each branch point is given by the minimum number of transpositions 
needed to build the conjugacy class corresponding to the partition^ For our case this formula is 
just a restatement of the Euler characteristic formula (|4ip . 

So each correlation function is a sum over Hurwitz numbers for maps to CP 1 with three branch 
points 

(59) 



tr(E„ Z^ k ) tr(S M Z® fc )) = (k\) 2 £ N c ^ Cov k (u, A, fj) 



There is a similar story for the Hermitian matrix model |24^56j. The relation with the complex 
matrix model discussed here is simple: just replace one of the profiles v with [2a] to account for the 
Hermitian matrix model propagator. There is then an identity that comes from setting t\. = 5k2, 
known in the literature \24\ [32] , 

\dZf NxN e -tr(S-4z5-^Zt)+E21l**tr(^)+tr(Zt») 
[dFf NxN e -trFirf + tr(,l-^SFt)<« (6Q) 

In the expansion of the Hermitian model in Feynman diagrams, /c-valent faces come with a coupling 
Sk while fc-valent vertices come with a coupling [24] . Graph duality in this model [5] just 
exchanges {t} for {s}. The identity with the Z model in the second line can be seen directly since 
the tr(Z^ 2 ) term allows each propagator to loop back to the holomorphic operator built out of 
Z's, which is then equivalent to the Hermitian matrix model operator in the Hermitian correlation 



9 There is a review of Hurwitz theory in [53]. The relation of these cut-and-join operators and Hurwitz theory to 
integrable hierarchies is summarised in |52| . 
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function. The difference in powers of S in the actions arises because there are now two edges in 
the Z model for every edge in the Hermitian model. The final line follows from the symmetry in 
the character expansion (|29p ; the identity between the first and third line was proved in |24j using 
the same techniques as used in Section [21^1 

Taking the limit described in Section [4.31 so that Y,\ — > £[2?] ~ pT^pl' w hich for the F model 
means restricting to only 4-valent vertices, the Hurwitz numbers become double Hurwitz numbers. 
Double Hurwitz numbers have two fixed branching profiles fi, v for branch points at and oo and a 
remaining arbitrary number of simple branch points (with profile [2] ) In this limit the partition 
function becomes the r-function for the Toda lattice hierarchy that appears in the Gromov-Witten 
theory of CP 1 [581 

co co k-2p/n \p 

Z({t},{t},{s})^Yl E E ' 1 k lJ tr K (A~T k ) tr(a„ (B-T k ) ^Sf 2] S,) 

fe=0 fj,,vhk p=0 

CO CO k — 2p/r) \p 

= E E E Sl J tr ^ ( A " 1 ) 0fc ) ^ ( B ~^ Cov ^ v >w>- Pi) 



V 

k=0 fi,p\-k p=0 



p times 



It should be possible to make this picture more precise: the Z and F models compute relative 
Gromov-Witten invariants for the topological A model on CP 1 with two points marked at and 
oo, cf. [6Uj . The holes in the F worldsheet corresponding to holomorphic Z operators wrap around 
while the holes corresponding to antiholomorphic Z^ operators wrap around oo. The parts of the 
cut-and-join operator T,\ then map to gravitational descendants T\.(w) of the Kahler class w. The 
full details of this correspondence require the technology of completed cycles [60] . 

The appearance of CP 1 \ {0, oo} is curious because of its appearance both as an auxilliary curve 
for the normal matrix model description of the c = 1 string |13| and in the topological B model 
set-up for the c = 1 string in [61] . The free energy of the c = 1 string at the self-dual radius is also 
known to agree with that for the topological A model on the conifold resolved by a sphere in the 
limit where the complexified Kahler class vanishes t — > [2]. 



5 The F model and the cell decomposition of A4 g ^ n 

The (Deligne-Mumford-compactified) moduli space M g , n of Riemann surfaces of genus g punctured 
n times can be extended to the space of 'ribbon' or 'fat' graphs by replacing each puncture with a 
boundary of length for i = 1, . . . n. This extended decorated moduli space M g , n x K+ has real 
dimension 6g — 6 + 3n. Due to the work of Mumford, Thurston, Strebel [62], Harer [63], Penner 
|35| and others, any point in this moduli space can be obtained by considering connected graphs 
with lengths assigned to each edge. 

10 Holomorphic maps onto CP 1 with just three branch points are special. Belyi's theorem [57] states that a non- 
singular Riemann surface is an algebraic curve defined over the algebraic numbers Q if and only if there is a holomor- 
phic map of the Riemann surface onto CP 1 with only three branch points. For the Hermitian matrix model, where 
all the ramification orders over one of the branch points are 2, the Belyi map is a special type called 'clean' |56] . 

11 Single Hurwitz numbers occur when there is no branching at the second point, i.e. v = [1 ]. The generating 
function for single Hurwitz numbers is obtained from the Z model with tk = Ski- 

12 The map to the variables in [55] is v., = tr(a M (A~ 1 ) 0fe ), p' v = tr(a„ (B~ 1 ) m ), q = si and /3 = 
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The Penner Hermitian matrix model |35] 

/ [dQ] H e tr«)+trlog(l-Q) = 




=2 k 



tr(Q k ) 



(61) 



gives a cell decomposition of M. g , n x Each graph of the Penner model with n faces and genus 5 
corresponds to one of the cells. The top-dimensional cells in Mg t n x are swept out by the lengths 
of the 6g — 6 + 3n edges of the Feynman graphs with only 3-valent vertices. Lower-dimensional cells 
in the cell decomposition arise when we shrink an edge, colliding two vertices into a higher-valency 
vertex. Because of the extra factor of —1 from the vertices with each lowering of dimension, the 
Penner model calculates the virtual Euler characteristic of Ai gn . The symmetry factors of the 
Feynman graphs account for the fact that M g , n is an orbifold space. 

Konstevich [3] adapted the 3-valent version of this model to give a generating function for 
the correlators of 2d topological gravity, which calculate intersection numbers on M g , n |H]- The 
couplings t n to the operators are encoded in a matrix Z by the transformation = ^tr(Z~ fc ). 
This constant matrix modifies the quadratic term in the Kontsevich matrix model 



In the expansion of this partition function, each Feynman graph with n faces and genus g can be 
written as an integral over the corresponding top-dimensional cell in Jvig^n x R™ using the Schwinger 
parameterisation of the propagators 



By integrating over the 6g — 6 + 3ra lengths p e of the edges of the graph the whole of the cell is 
covered. Separating out the integral over the boundary lengths R™ (which correspond to sums of 
the edges around each boundary) one is left with an integral over M. g>n corresponding to the closed 
string correlation function [3] (see [651 139] for concise summaries). 

The F-model only has vertices of even valency, which suggests that it localises on lower di- 
mensional cells in the complex for M g , n x R+- The maximal-dimensional cell for the F model, 
corresponding to all vertices of valency 4, has dimension 4g — 4 + 2n corresponding to the number 
of edges in the diagrams. This localisation on degenerate subspaces of the moduli space was noticed 
in [66} [67] when considering extremal 4-point functions in N = 4 SYM. 

Following the example of the Kontsevich-Penner Hermitian model of Chekhov and Makeenko 
[36] we will write each generic F model graph as an integral over a discrete version of a lower- 
dimensional cell in M g , n x R+ [ID]- in the continuum limit the Chekhov-Makeenko model reduces 
to the Kontsevich model. In the continuum limit of the F model we get a 4-valent model which is 
the complex analogue of the Kontsevich model. 

With discretisation parameter e rewrite A = \f~Ne eL and B = \/ r Ne eM , i.e. = \fNe eli and 
bj = VN e emj . The propagator (|27j) of the rescaled F model (|26p can then be written 



The sum is a discrete Schwinger parameterisation of the edge length for the propagator. Each 
summand comes from an edge of the F graph with integer length pe, corresponding to a different 




(62) 




(63) 




(64) 
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number p of bunched propagators in the dual Z model graph, cf. (|25p f^l The integer length t of 
the boundary of each face is the valency of the vertex of the dual Z model graph, or in other words 
the power of the operator tr(Z e ) or tr(Z' ). 

For an F model graph with V vertices, E edges and faces corresponding to Z vertices labelled 
by / (and colour index if) and Z' faces labelled by g (and colour index j g ) the contribution is 

cN v y rr 1 „ = cn v - e y TT -, 

^ 11 a it bj -N 11 el if enij , 

{ifUjg} hedges V •»« {if},{j g } E edges e Je 1 

£ oo 

{v}.{i 9 } r=1 P'-= 1 

c is the symmetry factor for the graph. The discrete sums for each of its E propagators give a sum 
over discrete points in an ^-dimensional cell of M g , n x WL. Each point corresponds to a different 
Z model graph with different numbers of bunched propagators between each set of vertices, given 
by the integers p r . A given Z model correlation function gets contributions from a finite number 
of graphs, so the closed string correlation function must localise on only a finite number of points 
in the moduli space. 

The propagator in its discrete Schwinger parametrisation (|64p has a continuum limit as e — > 

00 poo -I 

lim V £e -P^i+m,) = / dpe - p ( k+mj ) = ( 66 ) 

£ ^°fr{ Jo k + irij 

The model with this propagator arises in a double-scaling limit of the F model where we take N 
large ./V ~ Jj. Rescaling the F model matrix F = i/eG the action ()26|) becomes 

- ti(AFBF*) + N ZZi I tr(FFt) fc 

= -Netr \(1 + eL + 0{e 2 )) G (1 + eM + 0(e 2 )) G f j + Nix \eGG ] + e 2 \GG ] GG ] + 0(e 3 ) 
= -Ne 2 tr(G^LG) - Ne 2 tr(GMG f ) + Ne 2 \ tr(GG j GG j ) + NO{e 3 ) 

-> - tr(G* LG) - tr(GMG f ) + \ ti{GG ] GG ] ) (67) 
This G model is a complex analogue of the Kontsevich model 



/ 



[dGf NxN e -^LG)-tr(GMGt )+ ^G&G&) ( g g) 



It has the propagator identified in (|66p which is similar to that of the Kontsevich model 

Looking at an individual Feynman graph of the F model we can also see that only 4-valent 
graphs survive in this double-scaling limit. Suppose a graph has vertices of even valency 2k. 
There are V = vertices in total and E = | Y1T=2 ^kik = Y1T=2 edges. Setting N = 4j- 

and taking the limit e — > then the expression (|65p is only non-vanishing if there is a factor of e 
for each edge (cf. {MP ) 

iV^-s = e- 5 (70) 



13 Such a discrete metric on the moduli space also appears for the Hermitian matrix model in |68l 169] . where the 
integer lengths also correspond to the number of bunched propagators between the vertices. 
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This means that E = 2V, which is satisfied if only 12 is non-zero. In this case we get 



!To c ^ n 



{ifh{j a }r=l 



p r =l 



-ep r (h f +m,j g ) 



■ e n/„ + 

^ n l if +m,n 

{if},{j 9 } Z edges l i 



-Pr(k f +m jg ) 



(71) 
(72) 



This comes from the corresponding graph of the G model. 

The integrals over the worldsheet boundary lengths must be decoupled with care, since 
there is at least one relationship between the boundary lengths: the sum of the Z boundary lengths 
must equal the sum of the boundary lengths. For example in Figure [8] the three-point function 
has only two independent boundary lengths, not three. Once this is done, it is still not clear what 
quantity we are integrating over (a subspace of) M. g , n - 

The fact that the G model has only 4-valent vertices relates it to the BMN limit from Section I4~3] 
which arose from limiting the F model to 4-valent vertices. However it is still not clear what the G 
model is calculating in this context. The obvious link would be if in the continuum limiting process 
the discrete sum only got contributions when p ~ - = y/~N (BMN-length operators propagating 
between vertices of the Z model), but this is not the case. This issue is left for the future. 



6 Conclusions and future directions 

In this paper we have studied correlation functions that correspond both to tachyon scattering for 
the c = 1,R = 1 non-critical string and to a half-BPS sector in free 4d N = 4 super- Yang-Mills. 
In the Z complex matrix model the closed string insertions correspond to vertices of the Feynman 
diagrams. The Z model is precisely dual to another complex matrix model called the F model. In 
the F model the closed string insertions are now associated to the faces of its Feynman diagrams. 
This duality can be shown using character expansions, or by integrating in and out fields to see the 
graph duality dynamically, following the programme set out by Gopakumar [S]- 

Using the example of the Kontsevich model, the correlation functions of the F model can 
be written as sums over discrete points in subspaces of the moduli space of punctured Riemann 
surfaces. These discrete points correspond to ribbon graphs with integer-length edges. 

This complex matrix model duality could provide a prototype for understanding the AdS/CFT 
duality microscopically. It may be possible to rewrite (perhaps just free) M = 4 super Yang-Mills 
as a dual theory, where local operators and interaction vertices from J\f = 4 SYM correspond to 
faces of the dual Feynman graphs. The correlation functions of this dual theory would be easier 
to write as string moduli space integrals, following the Kontsevich schema. At non-zero coupling, 
summing over the interaction vertex holes should remove the D3 branes and alter the background 
to AdS 5 x S 5 . 

There is a long way to go to realise this goal. In contrast to the programme set out in [57 ] 155 ] 155] . 
here we have dropped the spacetime dependence of the N = 4 correlation functions to focus on the 
combinatorial index structure from which the non-planar expansion comes0 Really the sector of 
N = 4 we have studied just computes a metric on multi-trace half-BPS states and doesn't contain 
any spacetime information.^ More general correlation functions not only include more interesting 

14 Another way of introducing the spacetime dependence is discussed in [70] . 
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spacetime dependence but should also get contributions from the full moduli space of punctured 
Riemann surfaces, not just subspaces. 

As long as we keep the separation between holomorphic and anti-holomorphic operators, it 
should be straightforward to introduce other complex scalar fields into this duality! 15 ! The non- 
planar expansion of the free theory in terms of cut-and-join operators discussed in Section [3] follows 
through with little modification [21]. Allowing the scalars to be real, or introducing fermions and 
the gauge boson, would introduce more complication. 

From the string side, it is important to understand how the reduction to the c = 1 string from 
the full IIB string on AdS§ x S" 5 works. Both contain a Liouville direction and it is tempting to 
identify the R = 1 limit with the small radius limit of the bulk geometry corresponding to free SYM. 
The c = 1 string in this limit is known to have various topological coset descriptions [7TJ [72j [73] , s0 
perhaps it is possible to use the cohomological reduction techniques of [71] for sigma models with 
supersymmetric target spaces to reduce the full bulk coset. 

On the other hand, an alternative strategy would be to take a topological description of the 
c = 1,R = 1 string and try to include the full PU(2, 2|4) symmetry of free 4d M = 4 SYM. For 
example, it is known that the free energy of the c = 1, R = 1 string agrees with that for the t — > 
limit of the topological A model on the resolved conifold [2]. Once one understands how tachyon 
scattering is reproduced in that setting, one can think about how to include more of the spectrum 
of M = 4. An A model approach seems promising given that the correlation functions were shown 
in Section 14.41 to count holomorphic maps onto CP 1 . A tentative connection to the A model on 
CP 1 was made at the end of that Section. 

Focusing now on the F model there are several areas that merit further study: 

• A brane interpretation of the F-model would be welcome, perhaps along the lines of the 
relation between the Kontsevich model and the open string field theory of FZZT branes for 
2d topological gravity derived by Gaiotto and Rastelli [3]. Such an interpretation of the 
model has been discussed in [TSJ [76] . 

• What mechanism localises the integral over the moduli space to discrete points in Section [S]? 
Given that this is also a feature of the Hermitian matrix model |40| , does it always arise 
in free theories? 

• The Z model and the F model capture tachyon scattering for the c = 1 string at the self-dual 
radius, but they do not include all of the discrete states or the SU (2) symmetry at this par- 
ticular radius. Perhaps vortices appear in the F model as holomorphic and antiholomorphic 
operators tr(F k ) and tr(_F' fe ) like their appearance in the similar six-vertex model [77| . 

• The Toda integrable hierarchy structure of the c = 1 string at the self-dual radius has not 
been discussed here from the point of view of the F model. 

• An algebraic geometry interpretation of the F model might correspond to the limiting case 
discussed in [78] for the c = 1, R = 1 string. 

• The authors of [61] reproduced both the Kontsevich model and the c = 1,R = 1 tachyonic 
scattering matrix by considering non-compact branes in the topological B model in a deformed 

1 See Appendix [D] for a sketch of how this might work with two complex matrices. 
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conifold background. A model similar to the Z model was also studied in [79] for the c = 
1,12 = 1 string. Where does the F model fit into this picture? 



Acknowledgements: We thank Sasha Alexandrov, Robert de Mello Koch, Rajesh Gopaku- 
mar, Vladimir Kazakov, Sanjaye Ramgoolam and Volker Schomerus for discussions and comments 
on early drafts. We are also very grateful to the organisers of the Second Johannesburg Workshop 
on String Theory, during which many of these ideas were aired. 

A Examples for F model 

In this section we consider explicit examples of how the F model Feynman diagrams reproduce the 
Z model correlation functions. First we consider the F model expansion. 

A sine le 4-valent vertex (tv(FF^FF^)) can only contract with itself to form a planar graph 
that is dual to the planar 3-point function of the Z model, cf. Figure [8] and Section lA.il for the full 
expansion. 

For two 4-valent vertices (tv(FF^ FF^) tr(FF^ FF^)} there are several possible diagrams: 

1. Each vertex entirely contracted with itself is disconnected and corresponds to two discon- 
nected 3-point functions. 

2. Each vertex with just one self-contraction each corresponds to the lollipop planar 1 — > 3 
4-point function, see Section IA. 31 

3. The planar diagram with each propagator connecting the two vertices corresponds to the 
planar 2 — > 2 4-point function all connected in a loop, see Section [A. 41 

4. The torus diagram with each propagator connecting the two vertices corresponds to the torus 
2-point function, see Section [A. 21 

All the cases for a single 6-valent vertex (tr(FF^ FF^ FF^)} correspond to the collision of the 
two 4-valent vertices considered above: 

1. One way of contracting it in a planar way with 4 faces corresponds to the 1 — > 3 4-point 
function Y diagram, cf. Section [A.3i that arises from colliding the two vertices from case 2. 
above. 

2. A topologically distinct way of contracting it in a planar way corresponds to the 2 — > 2 
degenerate (tr(Z p ) ti(Z^ p+q ) tr(Z^ q+r ) tr(Z r )) from the collision of vertices in 3. above, cf. 
Section IA.41 

3. The vertex contracted with itself on the torus has two faces and is A4\ } 2, cf. Section IA.21 
This results from the collision of the vertices in 4. above. 

A.l M ,3 

The leading planar term of the three-point function in the Z model is 

(tr(Z tfc )tr(Z fcl )tr(Z fc2 )\ = k 1 k 2 kN k ~ 1 (73) 

\ / sphere 
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where k = k\ + fc 2 - This corresponds to the cut-and-join operator Er 2 i which splits tr(Z^ k ) into two 
pieces, cf. the analysis in [21 j . where there are two bunches of homotopic propagators, see Figure 
[8l The dual graph, drawn in Figure [8] has one four-valent vertex, two propagators and three faces 
(which correspond to the old vertices). Reading the term in the rescaled F matrix model (|26p 

|(tr(FFW)) = ^ £ (i^tg^tg 

= y E ((^^^0(^ Ft ") + (^ Ft ")(^ Ft -)) (74) 

U,*2 ,31,32 

Note that there are no non-planar terms in this correlator because of the configuration of the 
fieldsH Inserting the propagator (|27p and then Taylor expanding each one 

A li £ 2)Jl 5 j2 \ / 

_ AT ^ N k 1+ k 2 -2 ~ ^1+^-2 



2 „kl+k2l ) kiik2 fei k2ik 1 +k 2 

i,jl,32 fcl,fc2=l » il tl,i2)J fel,/C2=l W *2 i 

1 oo 

- [ir(vl- fcl - fc2 )tr( J B- fel )tr( J B- fc2 )+tr(>l- fcl )tr(yl- fe2 )ir( J B- fel - fc2 ) 



2 

fcl,fc 2 = l 



iV 



fel+fe 2 -l 



[tfcitfatfci+fa + ifciWfci^] fei*2(fci + /e 2 )iV fcl+fc2 ^ 1 + 4 E + t 2 k 1 t kl tk 1 } 2k\N 2k ^ 1 



2! 

k±<k 2 k\ 



This agrees with the expectation from (J73J) where we get contributions from the conjugate corre- 
lation function too. Note that the generating function splits into two pieces depending on whether 
k\ = fc 2 , in which case we get a factorial from the exponential in the Z action ([!]). 

A.2 Mi ;2 

The torus two-point function for the Z model is 



tr(^)tr(Z fc )) toms = (k\ (S [3] + E [2 , 2] ) \k)N^ = k [Q + Q 



iV fc " 2 (75) 



Here we've used the cut-and-join notation of [21]. The different cut-and-join operators correspond 
to bunching homotopic propagators into either 3 or 4 bunches, cf. Figure [9] for the two possibilities. 
The bunching of the propagators into three, the lefthand diagram in Figure [U yields a dual graph 
with a single six-valent vertex, three edges and two faces. Reading the appropriate term in the F 



6 In a Gaussian Hermitian matrix model (tr(M 4 )) does however receive non-planar contributions. 
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matrix model we compute the non-planar torus term for the six-valent vertex 
— (tr(FF^FF^FFU) =— V / Fj 1 F^ 1 F 1 - 2 F+f K 3 F^ 3 , 

3 \ V V torus 3 ^ \ n 12 n 13 33 11 /torus 

*1)*2,13)J1,J2,J3 

= ^ / P'i 7?ti2 \ / j^is pf ii \ / F i2 F ]h 

3 Z-^ \ ii »3 / \ J3 12 / \ h u 

11,12,13 Jl J2 J3 

at ,r ( y'-* .^m' tf 2 <tf 3 

3 11,^1,32,33 " " ^ " 

= T S ( ai 6 7 - - iV) 3 (76) 
Now Taylor expand, using the binomial theorem 



EA:Qt fc t fc iV fc - 2 (77) 



i,j k=3 x '13 k=3 

00 

rfc-2 



This agrees with the expectation from ([75 p . 

The bunching of the propagators into four, the righthand diagram in Figure El yields a dual 
graph with two four-valent vertices, four edges and two facei^l: 

(tr(FF^FF^) tr(FF^ FF^)\ 

2 • 2 • 2! \ / torus 

A/"2 , . 

8 ^ \ jl h j2 h j3 h j4 iz /torus 

«l,«2,»3i»4jl J2 J3J4 

There are two relevant non-planar contractions 

' pil ptJs\ / Z?«2 ptM / 77*3 ptM / 77*4 ptJ'A + / P»l /rtM / 77*2 Ptis\ / 77i3 /rtJ'A / 77*4 FtM / 7 o\ 
v Jl »4 / \ J2 13 / \ 33 *1 I \ 34 h / ~ \ 31 »3 / \ J2 «4 / \ J3 «2 / \ J4 U / V ' 

These both give the same contribution so we get 

V- 1 _ N*_ y y. (k - l\ iV^ 

4 2^,-^- 4 3 ja*&J 

00 



fc=4 ^ ' 



2 



(79) 



17 One might worry about disconnected non-planar graphs where each 4-vertex only contracts with itself. Fortu- 
nately this diagram is not possible, cf. the A1o,3 example. 
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A.3 M 0A 1 ->■ 3 

The leading planar term of the 1 — > 3 four-point function in the Z model is 

/tr(Z tfc )tr(Z fel )tr(Z fc2 )tr(Z fc2 )\ = kik 2 k 3 k(k - l)N k ~ 2 (80) 

\ / sphere 

where k = ki + k 2 + k%. This result is made up from contributions from the two cut-and-join 
operators Er 2j 2] an( i ^[3]j corresponding to the lollipop diagram with 4 edges and the Y diagram 
with 3 edges, as explained in |66} 12 lj . The lollipop comes from two 4- vertices of the F model with 
a self- contract ion each. The Y diagram comes from one of the planar ways of contracting a single 
6-vertex of the F model. 



A.4 M ,4 2^2 

For k\ + k 2 = k% + &4 and max{/cj} = k\ the planar connected 2 — > 2 four-point function is 

(tr(Zt fcl )tr(Zt fc2 ) tr(Z fc3 )tr(Z fc4 )\ = kik 2 k 3 k 4 {ki - l)N kl+k2 ~ 2 (81) 

\ / planar, connected 

For generic values of the {k{\ the Z graph has 4 bunched edges corresponding to Sr 2) 2] and the 
dual graph is the planar contraction of two 4- vertices of the F model with no self-contractions. In 
a degenerate case where k\ — k% = k 2 — k 4 there are only 3 edges in the Z model skeleton graph 
corresponding to S™ and the dual graph has a single 6-valent vertex. 

Note that this correlation function also receives contributions at order N kl+k2 and ]V fcl+fc2-2 if 
k\ = k% from disconnected two-point functions. 



A.5 M 



0,2 



The F model diagram dual to the planar two-point function requires special treatment. The dual 
graph (see the righthand side of Figure [8]) must be take using the plain Gaussian propagator ([24"]) 
with the quadratic interaction terms in the faces between the parallel propagators of the Z model. 
The F model diagram corresponds to taking these quadratic interaction terms daisy-chained with 
no self-contractions 



N 
~kJ 



tr(A~ l FB- l F ] ) 



k \ iV fc /-rrv- 1 1 




h ■ 3p 1p 

plain, daisy n ' : \„ =1 j ,• "%> 3p / . . 

plain, daisy 



„_-, ,■ u i P u j P x ' p iam 

Here ik+i = ii and similarly for jk+i- There are {k— 1)! completely equivalent ways of daisy-chaining 
the Wick contractions. Now inserting the plain propagator (|2T ~ 



x n e ^ - t e - M - tA ( tr(zl) tr(z,i) U re < 83 » 

p=i« P :ip p p l d 1 3 

In this final step we have used the formula for the planar two-point function. 

Note that because we are using the plain version of the F model there is no sum over topo- 
logically identical Z graphs with bunched propagators; here the Z correlation functions must be 
calculated separately for each k. 
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B Examples for C, D model 

In this section are calculations for the C, D model (|15p with quartic vertex and propagator (|16p . 



B.l M , 3 

For the specific 3-point function 

(tr(Z)tr(Z)tr(Z t2 )\ = 2N (84) 

\ / sphere 

exactly one C, D diagram contributes. This diagram is the same as (3) in Figure [6] except that 
there are only two quartic vertices. Proceeding with Einstein summation on ek, fk = N only 



2! E ( C zrn C ) 1 i 1 D hh D l 1 j 1 C e 2 i2 C ) 2 i 2 D hh D l 2 , 
h,i 2 ;ji ,32 

_ ^ei/i^nn ^/ie 2 ^jij2 ^eih^ji'h ^e 2 f 2 ^i 2 i 2 

2 . tt^ &71 tt^o 

\ E ^> = 2N ^ 



2 .4-'. Ojj a, 2 6f 

n,«2y J 



B.2 Mi, 2 

For the specific torus 2-point function 

(tr(Z 3 )tr(Z t3 )\ = 3iV (86) 

\ / torus 

the relevant torus C, D diagram has three quartic vertices 

31 E { Ce ^ C hh D fdi D l 1 j 1 C ^i2 C f 2i2 D f2j2 D l^^ ( 87 ) 

{ik},{jk} 

There are two ways of Wick contracting; one choice gives half of the result 

Jl E ( C ^ C L) ( C L C ™) { D hhD\ 2h ) (D\ xh D hj3 ) (C £212 C} 313 ) (D hh Dt 3h 

i,j 1 3 

C Relation of F model to Won model 



In this section we show that the Woo model \A2\ [33] can be directly related to the F model by a 
change of variables^ Take the model with v = —i/j, = N, so that the logM term is tuned 
away, and expand the exponentiated operators in the same way as in equation (|48p 



J[dM]^ N e- t ^+^^ t ^ M ^^ xdB-^xdMA- 1 ) (89) 



l(R)<N 



18 The proof in this section was carried out in collaboration with Hanna Gronqvist of the University of Helsinki. 
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The eigenvalues of M must be positive semi-definite for this integral to be well-defined. Shortly 
we will see how this condition is automatically implemented by the F model. Make the change of 
variables M = UDU^ where U is unitary and D is diagonal with eigenvalues mi, . . . > 

N 

Wjldmi A 2 (mi) e-*W £ X r(B^)xr(UDU^ A" 1 ) 

J i=l 1{R)<N 

= [ {[dm, A 2 (mi) £ Xn(B^) ^^f^ (90) 

i=l l(R)<N 

A(?rij) is the standard Vandermonde determinant. In the final line we have used the integral [22J 



[dU] U NxN X r(UXU^Y) = XR ^™P (91) 

We will now manipulate the F model in a similar way to get the same answer (|90p . The complex 
matrix F can be written with two unitary matrices U, W and a diagonal matrix D 



F = W\fDU^ F ] = UVDW^ (92) 

The eigenvalues mi, . . . m n > of the diagonal matrix D are the real, non-negative eigenvalues of 
FFK The measure is then 



m C e -tr(FFt)+iVEr= 1 j:tr[(A-iFB-iFty 

L J Yh X 7*1 

= [dU][dW][[drn i A /! (rni) e K 1 ^ k - x k L v > J (93) 

J i=l 

Character expanding the exponential with (|48p and using ()9ip on the [<iW] integral 

/n 
[dU][dW]Y[dmiA 2 (mi) e" tr{D) ^ X*(Ijv) XR^Wy/DU^ B^Uy/DW^) 

i=l l(R)<N 

= J [dU] n d m AHm) e~*™ £ dim - R XR{A ~ 1] ^R B ~ 1U) (94) 



i=l l(R)<N 

Next use (|91|) on the [dU] integral 



Jtk m<N dim n Rdna n R 

This already agrees with (|90p if we choose n = N. For the case n > N, compare (|95p with the 
character expansion of the F model (j29|) to see that 

„ n 

J H^AVi) e~ tv{D) XR{D) = [dim„,i?] 2 (96) 



Inserting this into (|9U|) we get agreement with 
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D Complex matrix model duality for two (or more) 

In this section we sketch how the duality might work for a V-type model with two N x N complex 
matrices X, Y, corresponding to two of the three complex scalars of free 4d N = 4 SYM: 



/ 



[dX] c [dYf e-^^-M^HEp^.M ^M 1 ,M 2 ,W } tr(aXMiy^)+ EMi ^ w l {MiiM2iW} tr(aXtM 1 ytM 2) 



The sum is over all the holomorphic single-trace operators built out of \x\ A°s and p2 Vs. a is a 
single /c-cycle a € [k] C S^ 1+ ^ 2 where k = \i\ + p.i- The trace with a permutation is defined by 

tr(a X^Y^) = X* 1 ■ ■ ■ X)^ Y- n+1 ■ ■ ■ Y^ 1+ " 2 (97) 

V ' Ml) ! a(ci) *Q(M1+1) Mmi+M 2 ) V ' 

It is unique up to conjugation a ~ p~ 1 ap for p £ 5 ftl x 5^ 2 so we only sum over conjagacy classes [a] 
for this relation. The couplings t, t can be encoded in a generalised Kontsevich-Miwa transformation 

*{Mi,l*2,[a]} = |Sym(a)n5 Ml xS M2 | tr ( Q ^C^) 

W,M} = IW^xg^l M" ^^ 2 ) (98) 

The matrices A,B,C,D do not commute and are not diagonalisable, unlike the single complex 
matrix case. For a single cycle Sym(a) = Some examples: 

*tr(X*) =i tr (^ fc ) *tr(tf»y») = tr(A 2 C 2 ) 

*tr(Y*) =itr(C fc ) ttr(xyxy) = £tr(ACMC) 

To get the dual model of -F-type the techniques of Section [2] using integration in-out look inappli- 
cable. Character expansions may work. A guess based on graph duality is 

J [dFf[dGf e" tr(FFt) " tr(GGt)+Efc ^w s ( fc -W} tr ( Q (afbf^i(afdg^(CGbf^3(cgdg^) ^ 

Each F propagator is transverse to an X propagator and similarly for G and Y . This guess has 
been checked for very simple two- and three-point functions. Here a is a single cycle permutation 
in Sk where k = J2i=i hi an d the coupling is defined 



{*<.[«]} ~ ISymMnn^feJ ( 10 °) 
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Figure 1: C,D integrated in. 




Figure 2: Z integrated out. 




Figure 3: F integrated in. 
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Figure 4: C, D integrated out; note that you can have any even- valency (FF^) P of F, F^ vertices. 
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Figure 5: Associations of A with dashed faces and B with dotted faces. On the left is the propagator 
and associated Oj,6j. On the right is an example six-valent vertex tr [(A~ 1 FB~ 1 F^)^, also with 
associated aj, bj. 




(5) o (6) 



Figure 6: (tr(Z 2 ) tr(Z) tr(Z"l" 3 )) graph duality, step by step. Graph (1) is the original Z model 
double-line diagram with three vertices and three propagators. In (2) C and D are integrated in, 
replacing the vertices of Z with faces of C and D. In (3) the Z propagators are shrunk to give the 
quartic vertices CC^DD^ of the C, D model. In (4) the quartic vertices are expanded in a different 
channel as propagators of F. (5) is the same graph as (4), just redrawn on the sphere so that the 
outer solid line of (4) becomes the inside central solid line of (5). In (6) C and D are integrated 
out; the faces of solid lines in (5) have been replaced by F/F^ vertices. Now faces of (6) bounded 
by dashed lines are associated with ^4's (corresponding to holomorphic vertices of (1)) while faces 
bounded by dotted lines are associated with B's (corresponding to antiholomorphic vertices of (1)). 
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Figure 7: Replacement of 2-valent vertices of the F model with plain ti(FF^) quadratic term 
(dashed lines), dual to faces bounded by parallel propagators in the Z model (thin solid lines), by 
the propagator of the proper F model (j23|) , which is dual to the parallel propagators of the Z model 
bunched into an edge of the skeleton graph. [All propagators are drawn in single-line notation here 
for ease of reading.] 




Figure 8: Z model three-point function and two-point function on the sphere with bunched prop- 
agators. The dual graph is drawn with a dashed line. 




Figure 9: The two different bunchings of propagators with no crossing for the Z model two-point 
function on the torus: three bunchings from Ew on the left and four bunchings from £[2,2] on the 
right. The dual graphs are drawn with dashed lines. The left dual graph can be considered as a 
degenerate case of the right graph when the two vertices of the dual graph on the right collide. 
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